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Abstract. Quasilinear perpendicular diffusion of charged particles in fluctuating electromagnetic fields is the focus 
of this paper. A general transport parameter for perpendicular diffusion is presented being valid for an arbitrary 
turbulence geometry and a plasma wave dispersion relation varying arbitrarily in wavevector. The new diffusion 
coefficient is evaluated in detail for slab turbulence geometry for two special cases: (1) Alfvenic turbulence and (2) 
dynamical magnetic turbulence. Furthermore, perpendicular diffusion in 2D geometry is considered for a purely 
dynamical magnetic turbulence. The derivations and numerical calculations presented here cast serious doubts 
on the applicability of quasilinear theory for perpendicular diffusion. Furthermore, they emphasize that nonlinear 
effects play a crucial role in the context of perpendicular diffusion. 
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1. Introduction and Motivation 

The influence of a turbulence on the spatial transport of charged particles plays a key role in a variety of heliospheric 
and astrophysical scenarios. The knowledge of the properties of the underlying turbulence is crucial for understanding 
the three-dimensional (anisotropic) diffusive particle transport in a collisionless, turbulent and magnetized plasma such 
as the solar wind or the interstellar medium. Of particular interest are the transport coefficients describing particle 
diffusion perpendicular to an ambient magnetic field. 

In spite of its long-standing importance in space and astrophysics, perpendicular diffusion has been an unsolved 
puzzle during many decades and a variety of studies have been carried out to achieve closure and to pin it down at 
a theoretical level. Models have been proposed using hard-sphere scattering in a magnetized plasma lOleeson llfltifll) 
and extended version of it developed on the basis of the Boltzmann equation 1,Tones [THflOII . but they are inapplicable 
to space plasmas where (electro)magnetic fluctuations trigger particle scattering. Other models are based on the field 
line random walk (FLRW) limit emerging, for slab turbulence geometry, from quasilinear theory (QLT; ,TokiDii [T9fif)l 
Jokipii & Parker flfltiflll . which has been considered and used in subsequent studies (e.g., Forman et al. IWTl Forman 
niTTTl Bieber & Matthaeus ITWIl . Although FLRW provides for a physically appealing picture, it has been shown 
recently that its applicability is questionable for particle transport in certain turbulence geometries, particularly 
those with at least one ignorable coordinate being the case for slab geometry (Jokipii et al. 1199Jt . Furthermore, 
numerical simulations, taking into account the magnetic nature of a turbulence, have shown that FLRW fails to 
explain perpendicular transport of low-energy particles (Giacalone & JokiDii 119991 Mace et al. I2()()()ll . nor is it clear 
that it actually accounts for the cross-field transport of charged particles. 

While perpendicular and parallel particle transport processes are necessarily distinct from one another in QLT, it 
has been argued for some time that parallel particle scattering can reduce perpendicular diffusion to subdiffusive levels 
if the turbulence reveals slab geometry only (see, e.g., Urch 119771 Kota & .TokiDii [?()()()! Qin et al. I2002b'l) . However, 
when the turbulent magnetic field has sufficient structure normal to the mean magnetic field, subdiffusion as seen in 
pure slab turbulence can be overcome and diffusion is recovered (see Qin et al. I2()()2;all . 

Recently, Matthaeus et al. proposed a promising model for perpendicular diffusion, also referred to as 

nonlinear guiding center (NLGC) model. The model takes into account different turbulence geometries and parallel 
diffusion while the particle gyrocenter follows magnetic field fines. Based on their model and its comparison with 
numerical simulations, it became clear that nonlinear effects play a crucial role for a more realistic understanding of 
perpendicular diffusion. 

Although more advanced models such as the NLGC approach have been developed, it is nevertheless instructive 
to provide a rigorous treatment of perpendicular diffusion in QLT. Related calculations by Shalchi & Schlickeiser 
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are valid for purely magnetic fields only. Furthermore, their calculations are restricted to specific wavenumber 
variations of the magnetic power spectrum and the turbulence correlation time. Detailed calculations do not exist 
allowing to consider quasilinear perpendicular diffusion in an electromagnetic turbulence. 

In view of this lack, it is desirable to derive a general QLT Fokker-Planck coefficient for an electromagnetic 
turbulence with arbitrary geometry. This and the detailed evaluation of the new coefficient for the limit of slab and 
2D geometry is the purpose of this paper. 

The structure of this paper is as follows: Section|2introduces the governing quasilinear equations of motions for test 
particles in fluctuating electromagnetic fields, and a general Fokker-Planck coefficient for a plasma wave turbulence 
with arbitrary geometry and arbitrary wave dispersion relation is derived. The limit of a slab turbulence geometry is 
considered in detail in section|3 There, perpendicular diffusion coefficients are presented for (1) an Alfvenic turbulence 
and (2) a purely magnetic turbulence. The limit of 2D turbulence geometry is considered in section 0] Numerical 
calculations and conclusions for both the slab and 2D geometry are presented in section 0 


2. Quasilinear Derivation of the General Diffusion Coefficient 

An often used standard approach for the evaluation of spatial diffusion coefficients is to calculate them from an 
ensemble of particle trajectories. To do so, the Taylor-Green-Kubo (TGK) formula is often applied. For the diffusion 
coefficient in x-direction, the TGL formula (KuboESZll reads 


Z 

= j di< Ua:(0)Ua;(C) > 


( 1 ) 


in the limit t —>■ oo, where Vx is the x-component of the particle velocity. The brackets < ... > denote an ensemble 
average over the relevant two-time distribution of particles. An analogous expression holds for the diffusion coefficient 
in y-direction. For a large coherence time the second-order velocity correlation function < Ua;(0)ua;(^) > must go to 
zero, and the integral in Eq. o approaches a constant value for t oo. Following the argumentation done by Kubo 
(see also Kota & Jokipii f^OOOII . equation o results from 


< Ax'^ >= 



I \ = 2 y - f) < Vx{0)vx{O > 


( 2 ) 


where < Ax'^ > is the mean square displacement of the particle position in time t. For t large compared to Eq. 
leads to < Aa;^ >= 2Kxxt, with the diffusion coefficient Kxx given by equation O- 

The use of the TGK formula is somewhat critical for situations where particle transport reveals rather an anomalous 
diffusion process than normal Markovian diffusion. In this case, the standard Fokker-Planck coefficient (anomalous 
transport law) 


< Ax'^ > 
2t 


(3) 


applies, where the mean square displacement scales more general as < Ax^ >oc . Depending on the exponent 7 , 
different diffusion processes can be taken into account: 7 = 1 for the diffusive regime, corresponding to a Gaussian 
random walk of particles; 7 = 2 in the superdiffusive regime, i.e. strictly scatter-free propagation of the particles, 
and 7 < 1 (e.g. 7 = 1/2, as shown by Qin et al. I2002bl and Kota & Jokipii [?000|l in the case of particle trapping 
(subdiffusion or compound diffusion). For 7 = 1 (normal Markovian diffusion), the anomalous transport law (|2I) is 
equivalent to the TGK formula o for the large t limit (see also Bieber & Matthaeus ITMTIi . 

In the context of QLT, the quasilinear perpendicular diffusion coefficients kxx and kyy can be written as 
(Schlickeiser 120021 Eqs. [12.3.25] and [12.3.26]) 


1 1 

fixx ~ 2 ; I^YY ~ 2 ^^.Dyy 


(4) 


where /i = U||/u is the pitch-angle of a particle having the velocity component U|| along the ordered magnetic field Bq. 
The subscripts X and Y denote guiding center coordinates, and Dxx and Dyy are Fokker-Planck coefficients (cf. 
Schlickeiser 120021 Eq. [12.1.29]) of the form 


Dxx = ^ d^< X{t)X*{t + 0 > ; Dyy = 3? / < Yit)Y*{t + ^ > 


0 


0 


(5) 
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They represent the interaction of a particle with electromagnetic fluctuations. The relations as given in equation 
correspond to the TGK formula o used earlier by Bieber & Matthaeus (jssai (see their Eq. [2]). In QLT, however, 
one has to perform an additional average with respect to fi. If the fluctuations are statistically homogeneous in space 
and time t, the velocity correlation functions in equation depend only on This can simply be taken into account 
in the Fokker-Planck coefficients by setting t = 0 , and the quasilinear relations then correspond, apart from 
the /i-integration, to the TGK formula cj. 

On the basis of Eq. 0 , QLT does not allow to consider subdiffusion or superdiffusion, since the derivations of the 
relations 0 are based on the assumption, that t is larger than the coherence time ^ (see Schlickeiser l2()()2l Eq. [12.1.17] 
and the comments following it). However, it was shown recently by Kota & Jokipii l|2flflfl|l that the Kubo formalism 
does not necessarily contradict compound (sub)diffusion if modihcations are applied. The consideration of anomalous 
particle diffusion (7 ^ 1), particularly of subdiffusion and an associated modification of the Kubo formalism, is beyond 
the purview of this paper and normal diffusion is simply assumed by applying the TGK approach. 


2.1. Equations of Motion 


The determination of Dxx and Dyy requires the knowledge of the equations of motion. According to Schlickeiser 
If^ll (see his Eqs. [12.1.9d] and [12.1.9e]), the perpendicular components of the fluctuating force fields can be written 
as 

gx=X{t) = -vcos(j){t)^/l - (SEy + (6) 

gy = Y{t) = -nsin(;ii(t)v'T^^7^-J - (SE^ - g-SBy) (7) 

dq tS^ \ c J 

where (f) denotes the gyrophase of the particle. Note that the Gartesian components of the fluctuating electromagnetic 
held, i.e. and j, y, are used and not the helical description, i.e. left- and right-hand polarized helds. 

For the further treatment of equations 0 and o, a standard perturbation method is applied. To do so, it is 
convenient to replace in the Fourier transform of the irregular electromagnetic held the true particle orbit x(t) by an 
unperturbed particle orbit, yielding 


(5B(: 


/ -l-oo „ 

* / d3/fc<5B(k,t)J„(lT)exp 


m['0 — + tkuvnt 


( 8 ) 


and an analogous expression for dE. The quantity J„(1T) is a Bessel function of the hrst kind and order n. The 
particle gyrophase for an unperturbed orbit is given by = (/>o ~ where the random variable (po denotes the 
initial gyrophase of the particle. Furthermore, the abbreviation W = k±RL\/l — is introduced, where Rl = v/fl is 
the Larmor radius. The relativistic gyrofrequency is given by H = qBg/{^mc) with m being the mass and q the charge 
of the particle, 7 is the Lorentz factor. The angle tp results from the wavenumber representation kx = k± cos ip and 
ky = k± sin Ip. With equation 0 , the equations of motion © and 0 can be manipulated to become 


X{t) = ^ Y I '^^^exp 


in [ip — (p(t)] -\- ik\\v\\t 


(9) 


X < - 




J„+i(W)e*’^+e-*’^J„_i(W) 


<5B|| -|— JnfW) i^Ey g—SBx'^ 


Yit) = J- E / '^^^exp 


n— — oo 


in [ip — (p{t)] + 


X < i 




J„+i(W)e*’^-e-*’^J„_i(W) 


SBu - 


^JniW) (SEx - I 


( 10 ) 


For the evaluation of equation 0 , it is convenient to consider now the nature of the electromagnetic turbulence. Here, 
the “wave viewpoint” is used by assuming that the turbulence can be represented by a superposition of N individual 
plasma wave modes, so that 

N N 

(5B(k, f) = Y^ (5B^ (k) exp(—zwjt) ; (5E(k, f) = Y^ (5E'^ (k) exp(—iwjf) 


(11) 
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Here, +irj(k) is a complex dispersion relation of wave mode j, where tOj^R is the real frequency of 

the mode. The imaginary part, rj(k) < 0, represents dissipation of turbulent energy due to plasma wave damping. 

Restricting the considerations to transverse fluctuations, i.e. • k = 0, and using Faraday’s law, the turbulent 
electric field can easily be expressed by the corresponding magnetic counterparts, yielding 


LOn 


[SB^-SBlky) ; 




[SBlk,-6Bm) 


( 12 ) 


Having expressed the electric by the magnetic field, it is also convenient to introduce now the Bessel function identities 

OiTl 

Jn-l{W) + Jn+l{W) = -^JniW) ; Jn-liW) - Jn+l{W) = 2j'^{W) (13) 

where the prime denotes the derivation with respect to W. With equations m and 0, the equations of motion 10 
and m can readily be rearranged, and one arrives at 


X{t) = Y I '^^fcexp 


j n— — oc 


in [ip — (pit)] + *(fc||'(;|| — LL!j)t 


(14) 


X Jn{W) 


a-^6Bl+b8Bl 

k± II 


- j;(w)| 


N 


j — 1 n— — oo 


(Tkexp 


in [ip — 4>it)] + ~ 


(15) 


X MW) 


a^SBl+bSBl 


+ iy^rMYsBiMw)^ 


where the following complex functions have been introduced: 




ujjk± 
V k^ 


b = 


-P 


(16) 


2.2. Velocity Correlation Functions 

Having determined the equations of motion, one can now proceed to calculate the second-order velocity correlation 
functions < X{t)X*{t -I- ^) > and < Y{t)Y*{t -I- ^) > entering equation 0. The procedure for the calculation is 
relatively lengthy, but can be carried out with simple algebra. For the sake of generality, the explicit time dependence 
of the correlation functions is included. At the end of this section, the limit < = 0 is considered to obtain < X(0)X*(^) > 
and < y(0)T*(^) >, required for statistically homogeneous conditions. 

The calculations for both velocity correlation functions are analogous, and the calculations are, therefore, restricted 
to < X{t)X*{t -I- ^) >. Multiplying equation itTHl with its conjugated leads to 


Xit)X*it + 0 = 152 


„ E E E 

U j n— — oo m—- 

x\MW)Jm{W) 


d^k / d^k exp(x) 


(17) 


kxkx 


aa 


(SBUBr) + bf ■ (SBiSB^’") 


:J*\ 


k±k± 


+bd*tx . {6Bi6Bi*) + ah* Y 


{5Bl5Bi*) 


kn 


+Mi- ^fmw)jL{w)^ 


-i^/TY?Jm{w)Mw) 


k± . 


k± 


i^-iSBlSB^ )+6-(<5H^dH|l ) 
fY.^SBlSBl*) + b*-iSBl 6 Bi*) 


kqi kii 


+(1 - M^MW)JLiW) ■ )| 

where the abbreviation 

X = i{n'ip — mip) — i{n — m)4>(t) + iW||(fc|| — k\\)t Y i(lO* — u!j)t — *(fc||'C|| ~^*j Y rnfl)^ 


(18) 
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has been introduced. The bar notation used over some quantities indicates that they have to be evaluated for wavevector 
k and time t + ^. A simplification of 113 can be achieved if one applies an average with respect to the random variable 
(jjo, the initial gyrophase of the particle. Using the relation 




27r 


d(j)o exp[^(n — m)(/>o] = Sn 


(19) 


where Kronecker’s 6nm = 0 for n ^ m and unity for n = m, and applying the ensemble average, equation 113 can be 
manipulated to become 


< 


oo oo 


Xit)X*{t + C)> = d^kexpix) 


j n— — oo m— — oo 


( 20 ) 


x<^ J„(IU)J„(1U) 


aa 


■* < SbIsb]* > +bb* < SBiSBi* > 


k±k± 


+ba*t^ < SBiSBi* > +ab*^ < SBUbI* > 

kj_ " k± II 


kx 

^k± 

kx 


+zv'1-mV„(1 U)J;(1U)^ 

k± 

-Vl-/iV„(W)j;(lU)-^ 

k± I k± 

+(1 - j;(iu) j;(W) < SBlSBf > 

where the quantity % is now given by 

X = — Ip) + ^v\l{k\\ — k\\)t + — ujj)t — i(fc||i’|| —oJ*+ nfl)^ 


a— < SBlSB\ >+b< SBiSB'^^ > 

< SBUBi* > +b* < SBUBi* > 


( 21 ) 


The next step consists of the assumption that Fourier components at different wave vectors are uncorrelated. 
Introducing the subscripts a and (3 for Cartesian coordinates, the ensemble averages of the magnetic field fluctu¬ 
ations can also be written as 


< SBiSB^p >=< SBi(k)5B^;{k') >= ^(k- k')Pi^(k) 


( 22 ) 


The uncorrelated state implies ip = ip^ W = W and coj = ujj, and the velocity correlation function, equation 1201 ), 
reduces to 


z p 

< X(t)X*(t + c) > = X/ / e'Kp[2Tjt - t{k\^v\\ - u;j^R + nil)^ + Tj^] 

x|j^(VF) 


(23) 


h 

an* — -I- bb* _: 

_L 




+*V'W^J„(IU) 4 (IU)^ 




Ik 


+(l-M")^[ 4 (IU)]"p,j 


Analogous to the calculations presented above, one obtains for the velocity correlation function in K-direction the 
expression 


2 cxi p 

< + ^) > = ^ / d^k exp[2Tjt — i{k\\v^\ — 


(24) 


j n— — cx 


i , 

ki. k^. 

- 

< JniW) 

aa*^P||ll+66*P4 + ^ 

ab-Pi^+lm-p’P 
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The advantage of having applied the Bessel function identities m to the equations of motion d and d is now 
obvious. Both velocity correlation functions and, therefore, nxx and kyy are expressed as a sum of three terms. Each 
contribution includes either J^{W), Jn{W)J^{W) or [JniW)Y■ Each term is accompanied by a specihc factor which 
contains the components of the magnetic correlation tensor and, therefore, the complete information about 

the turbulence geometry. 

Furthermore, note that the velocity correlation functions reveal an explicit dependence on time t. This is a con¬ 
sequence of the dissipation of turbulent energy due to plasma wave damping. To demonstrate this, equation d 
is considered again but, for simplicity, without the concept of a superposition of different wave modes, so that 
(5B(k, t) = (iB(k) exp(—iwt) where uj = ujr -I- zT with the dissipation rate T < 0. The corresponding turbulent 
energy can then be expressed as 

< SB{'k,t)SB*{'k,t + ^) >=P(k,t,e) =< SB{'k)6B*{'k) > =P(k,C)e2^* (25) 

yielding, for a dissipative system, the well known relation 

(26) 

where 2r is the rate of energy dissipation due to collisionless plasma wave damping. According to the TGK formula 
CJ, one has to evaluate the velocity correlation functions for the condition t = 0, e.g. < A(0)A*(^) >, since the 
fluctuations are assumed to be homogeneous in space and time, and the variation in t vanishes. Hence, the calculations 
presented below do not take into account the contribution 2Tjt. However, such a contribution may be of interest for a 
more general treatment, whether in QLT or in an extendet version of the NLGC model, which is not considered here. 


2.3. Fokker-Planck Coefficients 

Having determined the velocity correlation functions in the previous section, one can now proceed and evaluate the 
Fokker-Planck coefficients as given by the relations in equation with t = 0 in the velocity correlation functions. 
Upon substituting the expressions ra and (1^ . one obtains 




j n— — OQ 


J^(lU)Fi, +ZJ„(1U)4(1U)G^^ + [4(lU)]'ff 


(27) 


and an analogous expression for Dyy, but including different auxiliary functions Fy, Gy and Hy. The latter and 
their corresponding counterparts for Dxx read as follows: 


. h V 

p3 — — J^hh*p3 _|_ ^ nh*P^ +hn*P^ 

ry — aa ^2 U||| + 00 -f- yao + oa 


(28) 


K = 1-^11 + bb-n, + A [<.!.•< + bo.-Pi 


v\\ 


(29) 






(30) 


Gi = -\/rr 




i-'‘')inpbPi,-b-Pi 


\\y 


^x 




(31) 

(32) 


p3 

Hy 


(1 _ !^p 3 

p ) ki n 


(33) 















O. Stawicki: On quasilinear perpendicular diffusion 


7 


The integration with respect to ^ leads to the complex resonance function, 

OO 

TZj — J [“*(^ 11^11 ~ 


Tj + z(fc||U|| — ujj^R + nfl) 

+ ih'^w - ^j,R + 


(34) 


which describes interactions of the particles with the plasma wave turbulence (remember that Tj < 0). 

The coefficients Dxx and Dyv, as given by equation can also be summarized to yield a net coefficient, i.e. 
D± = Dxx + Dyy. According to equation 0, one can then define a net perpendicular diffusion coefficient 


+1 

Ky = KXX + KYY = 2 / 


(35) 


with 




3 


Jl{W)Fi+iUW)J'^{W)&^ + [J’^iwfHi 


(36) 


The auxiliary functions G(l and are then a simple superposition of the corresponding functions given by 
equations ra to iISSIi, resulting in 

fl 1 t-iT 1 fl T 1 (37) 


Fi = «a*P|f||+56*(Pi, + P4) + — ^ 


h 4- Jc 


L 


Jc P^ —I— Jc P^ 
-I- >^y^y\\ 


Gi = 


Vl - F 


2 r 




pj 

C-Ty 


.)-■( 






(38) 


Pi = (1 - ^2) P|f|| (39) 

The coefficient (P|l. one of the main results of this paper and presented here in this general form for the first time, 
enables one to calculate the perpendicular diffusion coefficient 13511 for an arbitrary turbulence geometry. Moreover, 
it allows to evaluate kj_ for a turbulence consisting of transverse wave modes with dispersion relations depending 
arbitrarily on wavevector. 

Whether for numerical or analytical calculations, a further treatment of its auxiliary functions to 13911 requires a 
certain representation for P;(^. In view of the structures of equations 13711 to 13911 . it is clear that different representations 

for P^^ will alter the underlying mathematical and physical structure of D± and, therefore, k±. 

Here, a representation is chosen commonly used in the literature. Following, e.g., Lerche & Schlickeiser 1200111. the 
nine components of P^^ can be expressed as 


Pi^iky,^)=Ayky,^) 




kakp j. , ki, 

-j^+ia^{ky,kii)ea,0,— 


(40) 


where the real quantity denotes the magnetic helicity, Sa/s is Kronecker’s delta and Ca/ji/ is the Levi-Civita tensor, 
and is the wave power spectrum. With the representation igni) for one can now proceed with the evaluation of 
equation For this, the appropriate components of (P|) are substituted into equations and (TTHll . Making 

use of equation , one arrives at 


Ei 

A-5 


iak±-bki\y ^ 


1 - 


luil 




2 /rA:|| 


(41) 


A^ 


k 


b- 


V k^ 


1 - 


^j,R 


(42) 
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where a and b are now real functions, i.e. 




n ’ vk"^ 


(44) 


with = ojjLj*. The magnetic helicity cr-i enters the coefficient D± only via G^j_, which reveals a purely imaginary 
character, whereas and are real helds. An evaluation of equation (IS51) involves, therefore, only the real part of 
the resonance function m, i-e. 


= - „ 
J t^2 


+ (^PII - ^j,R + 


which is a positive-definite entry in the net diffusion coefficient k_l, since T^ < 0. 


(45) 


3. Diffusion Coefficients for Slab Geometry 

Although is has been found that the solar wind turbulence is often dominated by its 2-D modes and has only a small 
fraction (say ~ 20%) of its energy in the slab contribution (e.g., Matthaeus et al. 119901 Bieber et al. I1994II . it is 
nevertheless instructive and desirable to have a solid treatment of quasilinear perpendicular diffusion in pure slab 
geometry. For slab geometry, the wave power spectrum can be given by 

A^=5^(fc||)^ (46) 

To examine Hj. for the slab geometry limit, the following approximation process is applied to equation (TTHll : the 
argument W = k±v±\/l — fj? of all Bessel functions is assumed to be much less than unity, since k± 0. Then, 
/c_L is small compared to fey, but not equal to zero. This leads to a “quasi”-slab model. Finally, the limit /cj_ = 0 is 
considered. 


3.1. Fokker-Planck Coefficient for Slab Geometry 

To start with the approximation, it is convenient to rewrite in equation ll^ the n-summation. Introducing, for 
illustrative purposes, an arbitrary function I{n), one gets 


I{n) = 7(0) -b ^[/(-n) -b /(n)] = Y 


(47) 


r—dzl n—0 


n—0 r=±l 


where the last step indicates that the sequence of summation may not be switched. Based on equation lITTIl . the 
coefficient can be written as 


= wYl Z! / d?knyrn) 


(48) 


0 j—:tlr—:kln—0'^ 

J^{W)Fi{rn) + i.UW)JUW)Gi{rn) + [j;(lF)] "i7i(rn) 


The argument occurring at the resonance function and the auxiliary functions indicates that one has to change in the 
corresponding equations the quantity n to rn. The Bessel functions are not affected by the new sum index r. For slab 
geometry, the asymptotic 

is used, representing the Bessel function for small arguments W —> 0, with r(a;) being the Gamma function. 
Substitution of expression 133) into the coefficient 133 results in 


= EE 

0 j=±l r=±l n=0 ‘ 


d^k 


TZj (rn) fW 


r 2 (l -bn) 


n (rn) + iG^j_ (rn) 


W 




(50) 
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A further treatment of equation CT requires a closer inspection of equations lUTl) . (Ba and (gSl- They can be 
represented as polynomials of different orders in the ratio njW'. the field F]_ is a second-order polynomial in njW^ 
whereas and are polynomials of zeroth-order. The term in the brackets of equation (PI then yields 


F]_ (rn) 


tG^j_ (rn) 


W 




A^(k±,kii) 



(51) 


where the coefficients ai, 02 and 03 read 
ai = 2(1 -m')|a 


02 


= \/2ai 


(fc||U|| -2ujj^R+\ujj\)r ^ / ujj,Rk^ 


vk 


vk'^ 


(52) 

(53) 


03 = 


i\k\\ 


vkP- 






vk"^ 

With equations (ICTi to (IMIl . the perpendicular diffusion coefficient (IKUll can be rearranged to become D±_ = £>±,1 
+ -D±, 3 , where 

2n-2 


(54) 




4£2 

U j=±lr=±l' 


n—0 


D±,2 = 


dk^a2A^k±,ki\)'^', 


Tij{rn) fW 


9 r 2 

0 j^±l r^±l 


V 


n—0 


^ E E / rffcW^(fc±,fcii)E 

0 j=±l r=±l d 


T^{l + n) V 2 
TZjirn) fW 


^ r 2 (l + n) V 2 

n—0 ^ 


(55) 

(56) 

(57) 


For the limit of pure slab geometry, the n-summations in the integrals are subjected to a closer inspection. The n = 0 
contribution vanishes in any case for and Dj_ 2 - For W oc fc_L = 0, it is clear that only n = 1 can contribute to 

the sum appearing in Dj_ i. Concerning Hq 2 , all contributions vanish, and one obtains 


V 2 T^jirn) fw\ 
L THl + n)[2) 


00 


; E 

n—0 


TZ-i (rn) 

71 _ — _-_ 

r2(l + n) 


The evaluation of the sum in equation leads to 


^^r 2 (l + n) V 2 J 


TZj{n 


0 ) 



(58) 


(59) 


The second relation in equation lEHli implies D ±^2 = 0. Consequently, the magnetic helicity does not affect QLT 
perpendicular diffusion of particles in a slab plasma wave turbulence. Substitution of the left-handed expression of 
equation (EHI into (CT . and making use of ( 1 ^ and the wave power spectrum (gni, yields 


D 


slab 

±A 


= ^a-M^)E E 

0 


dk\\g^{k\\)TZj{r) / dk± 


S{k±)k'i 

kf,+k^, 


= 0 


(60) 


Apparently, i 4 j _,3 is the only nonvanishing contribution for slab geometry. Rearranging some terms in 03 , equation 
(PI , one arrives at 


D 


slab 

± 


j~)slab 

-^±,3 


2'nv^ 

R 2 

^0 


00 

■ E / dk\\g^{k^\)nj{n = Q) 

f=±l-oo 

{k\\v^^-^3,R)'^ r|(fc||) u 

v'^k?, r;2/j2 


{\ujj \ -UJj^r) 


(61) 


Equation m is the general FPC for perpendicular diffusion of particles in a slab turbulence consisting of transverse, 
damped plasma waves. The wave power spectrum g^, the real frequency uJj^R of the plasma wave mode and the 
corresponding dissipation rate Fj depend arbitrarily on fc||. Equation m allows to derive a general diffusion coefficient 
f^siab only for the plasma wave viewpoint, but also for the so-called dynamical magnetic turbulence approach. 
Perpendicular diffusion coefficients for slab geometry and for these two approaches are derived in the next two sections 
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3.2. Plasma Wave Turbulence Approach 

After having considered the slab limit of the general Fokker-Planck coefficient (|i^ in some detail, the corresponding 
diffusion coefficient is derived in this section. For this, the pitch-angle integration as given by equation is 
performed and then, if necessary, the integration in wavenumber space. Substitution of equation inn into results 
in 


K 


slab _ 

± — 


1 

2 


+1 

J 


-1 



dfclig^F, [/i +h+ h] 


(62) 


where the functions /i, I 2 and I 3 are pitch-angle integrals. Furthermore, the resonance function 123) has been used. 
The pitch-angle integrals can be solved analytically, yielding 


h 


+s_ 

J 

"-B+ 




arctan 



-I- arctan 



h 


v^ku 


+ B_ 


dx 


' -B+ 


q + X^ 



arctan 



+ arctan 



(63) 


(64) 


h 


r X+UJj^R 

v^kf. J r2 + 

" -B+ 


1 


-k r|) - ln{Bl + F^) 



arctan 



+ arctan 



(65) 


where the integration boundaries are given by B± = fc||U ± Concerning 

t^j,R (1 + P‘j/^j,R) — 1 in equation lIHKll . it is obvious that I 3 becomes negligible for 


the term \Ljjj\ — l^j^r = 
the condition |Fj| «C |wj,_r|. 


This is an often used assumption within quasilinear theory to derive the dissipation rates of plasma wave modes, and 
the approximation I 3 = 0 is considered as reasonable for further progress. 

With equations (P|l and the quasilinear perpendicular diffusion coefficient for a plasma wave turbulence with 
slab geometry reads 


K 


slab 

± 




arctan 



-k arctan 



( 66 ) 


where it was assumed that < 7 ^(fey) = g^dfcyl), rj(fc||) = rj(|A:|||) and ujj^R{k^^) = ujj ji{\k\^\). Furthermore, the fact was 
used that arctan(a:) is an odd function in x. As a consequence of this, the dissipation rate Fj is now given as a positive 
quantity, meaning that Fj > 0 in equation ra . 

Further evaluation of the wavenumber integral in equation II 66 II requires specific wavenumber variations of gT tOj^R 
and Fj. Especially the dissipation rate Fj can usually be given as a quite complicated function not only in fcy, but 
also in the so-called plasma fip.e = Sirnp^ekRTp^e/B q. In the latter, the subscripts refer to protons and electrons as 
plasma constituents having the corresponding temperatures Tp_e and the number densities rip^e- The temperatures of 
the plasma electrons and protons can differ quite significantly, giving rise to instabilities (e.g., Gary initial) . It is well 
known that such a plasma configuration with different temperatures of the plasma components can then result in a 
variety of wave modes with substantially different tOj^R and dissipation rates F^. 

For instance, consider the Alfven wave mode with real frequency tOj^R = jvAk\\. A typical representation for the 
damping rate of this wave mode can be derived by using a quasilinear Vlasov theory code (e.g., Garv IlihJHIl . Doing 
this, one obtains 


r,(fc||) 


0.60/3°'36Op 



1 . 54 / 30‘’3 

exp 



(67) 


where cOp and flp are the plasma frequency and the plasma proton gyrofrequency, respectively (Stawicki et al. 1200111 . In 
view of O) . it is clear that the wavenumber integration in equation has to be solved numerically. This is beyond 
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the scope of this paper, and, for illustrative purposes, a simplified version of the damping rate O) is demanded. For 
this, the exponential expression in Eq. is roughly neglected. Furthermore, it is assumed that /3p is such that the 
exponent of the power law is equal to unity. The two crude simplifications lead to 


ri(fc||) =rj.or^Afc|| (68) 

where Fj^o is a constant. Substitution into leads to 




slab 

± 


4,'Kvrj 

R2 

^0 



arctan 


^ 1 + vaIv \^ 


+ arctan 


(1^) 



(69) 


where it was assumed that r_|_ g = r_ g = Fg and = g = g- Furthermore, the quantity rj = TqVa/v has been 
introduced. Apparently, the expression in the curly brackets of iinni) can be shifted in front of the wavenumber integral. 
It has to be stressed that this is a consequence of the crude approximations made to simplify the dissipation rate 
Fj is linear in fey and, therefore, wipes out the wavenumber dependence of 

Taking into account that va/v 1 and neglecting the corresponding terms in equation (I69II . one obtains 


CXD 

^siab ^ ^2-77 arctan J dfc|| (70) 

A closer inspection of equation rni) . and its general version, equation ra . results in the following important obser¬ 
vation: vanishes for a dissipationless plasma wave turbulence, i.e. Fj = 0. Consequently, for a static turbulence, 

charged particles can not move in perpendicular direction and are trapped to a line of force of the background magnetic 
field. 


3.2.1. Dissipationless Alvenic Turbulence 

Tlj{n = 0) = TrS{k\\vg.-jvAk\\) (71) 

where the Alfven dispersion relations ujj^r = jvAk\\ was used. For F^ = 0, the second and third term in the brackets 
of equation m vanish in any case. The first term remains as the only one for the following two derivations. 
Derivation I: In order to obtain the wavenumber integration is first performed, and then the /r-integration. 

Making use of equation GB, one obtains for the FPC m the expression 

OO 

^0 .=±i_i ^0 i=±i 

The argument of the wave power spectrum g^ indicates that is has to be evaluated for A:|| = 0. Equation 18511 then 
yields 


+1 


^slab 


= - d^D 


slab 


= 


(o) + <?-(o)] 


^0 


dg, 


\g-jvA/v\ + \g + jvA/i 


(73) 


The integrand can be approximated by 2vaIv and 2g for the intervals 0 < g < va/v and va/v < p. < 1, respectively, 
and one arrives at 




(74) 


In the last step of equation iTTl) it was made use of g'''(0) + g (0) = SB‘^{0)L/4tt, with L being a scale length 
providing for the right normalization. The neglect oi v\/v^ 1 corresponds to a vanishing irregular electric field and, 

therefore, results in a purely magnetic turbulence. Equation m agrees with earlier quasilinear results interpreted as 
non-resonant field line random walk in a magnetostatic turbulence ('.Iokipii fl966l Forman I1974II . 

Derivation II: The second derivation is based on the fact, that one can easily switch the order of integration: first, the 
/i-integration is performed, and then, if necessary, the integration over fc||. Substitution of the resonance function GB 
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into Eq. results in 
2 ^ 

^ X! / dkiig^{ki\)ku^6mv^i-jvAki\)mvfi-jvAk\\f 
0 .=±1 J 

— OO 
OO 

/ dfc||5^(/c|i)fc||'M(Ai - juA/w)(Ai - 

X_LI 'J 


Stt^u 

R2 

^0 ^=±1; 


Applying now the pitch-angle integration, according to Eq. ((23), one obtains 

- 1-1 


^slab 


R2 

^0 ^=±1 


^ j J dfi5{g - jVA/v){ti - JVA/vf = 0 


( 75 ) 


(76) 


Therefore, the perpendicular diffusion coefficient vanishes for a dissipationless turbulence. As expected, the result GSl 
then agrees with the general diffusion coefficient ra . 


3.3. Dynamical Magnetic Turbulence Approach 

In this section, derivations of quasilinear perpendicular diffusion coefficients for the so-called dynamical magnetic tur¬ 
bulence (DMT) model, introduced by Bieber et al. II19H4II . are presented. Within the context of DMT, fluctuations are 
assumed to be purely magnetic. To take into account the dynamical behavior of such purely magnetic fluctuations, 
Bieber et al. llT!Mll defined two models: the damping as well as the random sweeping model. It is relatively straight¬ 
forward to derive on the basis of the more general FPC 16111 for these two models, and both are considered in 
turn. The limit of vanishing turbulent electric fields can easily be achieved by setting = 0 and Tj = 0 in Eq. 
EH), initially derived for the plasma wave viewpoint. The only modification concerns the real part of the resonance 
function m- Furthermore, since the concept of a superposition of individual wave modes does not apply anymore, 
the j-nomenclature is dropped. 


3.3.1. Damping Model 

For the damping model, Bieber et al. dnna suggested a dynamical behavior of the turbulent energy being of the form 


< SB{k,t)SB*{'k,t + ^) >= P(k) exp(—t'cO (77) 

and the resonance function m has to be modified in this respect. What has to be done is to set Tj = 0 in the upper 
part of Eq. m, and then multiply the exponential expression with exp(—The contributions resulting from the 
unperturbed particle orbit still hold. Having in mind that one has to take the real part of equation (EH (see Eq. iTOll 
and the comments following it), the resonance function for the damping model reads 


OO 

3^7?. — d^exp [—*(fcp|| + nil)^ — = 


+nny 


Following Bieber et al. ra, the rate for turbulent decorrelation, Uc, is assumed to be 


(78) 


Uc = amVAk\\ (79) 

where the parameter 0 < am < I allows adjustment of the strength of the dynamical effects. The case am = 0 then 
represents the magnetostatic limit, am = 1 describes a strongly dynamical magnetic turbulence. With the decorrelation 
rate (PI and the resonance function one readily obtains for Eq. EH the expression 


= J dfc||g(fc||)M^7^(n = 0) 

— OO 


SttuC p3 
Bi 



(80) 


where C, = amVAlv. Furthermore, g{ky = ( 7 (|A:|||). With equation (1851) . one obtains 


^slab _ 

Kj_ — 



-1 


SttvC 


dki 


gih) 


dp 




1 -I- p'^C 


0 


0 


(81) 
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where the /^-integration can be carried out analytically, yielding 


^slab 


8ttv( 

L 


1 — C arctan 


dki 


5 (^ 11 ) 


(82) 


Obviously, for the limit am = 0 (magnetostatic limit), vanishes, implying that the particle remains tied to the 
background magnetic field. Shalchi & Schlieiser (EEl came to the same conclusion. 

3.3.2. Random Sweeping Model 

For the random sweeping model, Bieber et al. used a Gaussian dependence for the turbulence decay, i.e. 

< dB{k,t)dB*{k,t + 0 >= P{k)exp{-i^^^^) (83) 

This changes the governing resonance function m significantly, and one obtains 

(A:||U|| -|-nO„) 


diTZ = j d^exp [—i(fc||n|| + nQ)^ — exp - 




(84) 


Making use of equation EJ and setting uij = Tj = 0 , one arrives at 

1 


1 


+1 


= 2 / 


f ,, g(fc||) 


dk\\ 


0 


={fc||) 


d/i/r^ exp — 


4z^c(A|i) 


-1 00 

The pitch-angle integral can be solved again analytically, resulting in 


slab _ 8^ 2 /. ^ 

d2„. / 7,2 


Blv 


kl 


0 


V^Erf 


k\\v 


k\\v 


2vc{k\\)) Vc{k\\) 


exp 


Y 

4i^c(*ll) / 


(85) 


( 86 ) 


where Erf(a;) is the Error function. So far, equation (I 86 II is the general representation of the perpendicular diffusion 
coefficient in a slab turbulence obeying the random sweeping model. Upon using for the decorrelation rate the same 
dependence on fey as for the damping model, i.e. Vc = amVAk\\, one obtains 


^slab 




0FErf(^)-^exp(-^ 


dk\ig{ki\) 


(87) 


Again, vanishes for the magnetostatic limit am —^ 0. 

4. DifFusion Coefficient for 2D Geometry 

The previous section offered detailed insight into the slab limit of the general EPC and its associated transport 
parameter Perpendicular diffusion coefficients are presented for the slab limit for both the plasma wave view¬ 
point and the dynamical magnetic turbulence approach. In this section, the evaluation of the perpendicular diffusion 
coefficient for the 2D contribution of the turbulence is presented. Eor simplicity, the calculations are restricted to the 
damping model of the DMT description (see also Section fOll . and the more general case of a plasma wave turbulence 
is left as an exercise to the interested reader. The calculations presented in this section generalize the derivation 
presented by Shalchi & Schlickeiser II2()()4II . since their approach is restricted to a simplified magnetic power spectrum 
and turbulence decorrelation rate. 

For 2D turbulence, wavevectors are perpendicular to the mean magnetic field, and the wave power spectrum can 
be given by 


A^{k±,ki\) = g{k±) 


sih) 

k± 


( 88 ) 


Upon substituting the expression lIHHll into equation and making use of the Bessel function identities the 
FPC for 2D geometry can be cast into the form 


= 


2ttv 


Bl 


2 00 




dk± 


0 


g{kA.)vc{ki.) r (1 - /r^) 

v'likx) + vP'VB 2 




(89) 
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Fig. 1. Numerical solutions of Eq. Ml representing the behavior of /(C, z) for three different values of z 


where cr = 0 was assumed. The corresponding resonance function has already been inserted with Vcik±) being the 
decorrelation rate in normal direction (cf. Eq. (I78II 1. According to equation 18511 . one has to perform the /r-integration 
to obtain the corresponding diffusion coefficient . The integration can still be carried out analytically, and the 
detailed calculations are presented in Appendix El There, it is shown that the diffusion coefficient for 2D geometry 
reads 

OO 

= J dk±g{kj_)I{C,z) (90) 

^ 0 

with 





cosh{29z) 
sinh(z7r) cos 9 


(1 — cos(26*)) sin(2^ cos 9) 


2( cos 9 cos(2^ cos 9) 


(91) 


where the abbreviations ^ = k^R^ and z = r'c/fl are introduced. Furthermore, Rl = v/fl is the Larmor radius. 
Equation (EOll is valid for a power spectrum and a decorrelation rate varying arbitrarily in wavenumber k±. The 
integral representation results from the /i-integration and has to be evaluated for further progress. Unfortunately, 
an analytical solution for this integral does not exist, and any progress requires numerical treatment. In order to obtain 
some insight into the behavior of the function /(C, z), Figure ^ shows numerical solutions of Eq. 19111 as function of C 
for three different values of z. For illustrative purposes, z is assumed to be constant in k±. 

The limit C ^ 1, implying Rl <C leads to an instructive, analytical solution for equation 19111 . To show this, 
small argument approximations for the circular functions are used, i.e. 


cos(2^cos6*) ~ 1; sin(2^cos6>) ~ 2^cos0 (92) 

and inserted into I{^,z). Partial integration then results in 

/(C«l,z) = -^ (93) 

Consequently, one obtains 

OO 

(94) 

0 

for the limit Rl ^ where the relation has been used. An eyecatching feature of equation lIMI) is the 

term including the dimensionless product Tcfl. It is formally the same as those derived by Forman et al. 1197411 for 
field line random walk in a slab geometry and, more recently, by Bieber & Matthaeus CM3. 

For the magnetostatic limit, 0 or z = 0, it can be shown easily that zt^^ ^ oo, by using small argument 

approximations for the hyperbolic functions. However, for slab geometry, the diffusion coefficients approach zero for 
z^c ^ 0 (see section EJ. 
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Kinetic Energy [MeV] 

Fig. 2. Numerical solutions of equation and representing the diffusion coefficients and , respectively, 
as functions of the particle kinetic energy. The calculations were performed for protons (p) and electrons (e“). 

5. Numerical Calculations and Conclusions 

To demonstrate the potential and flexibility provided with the new perpendicular diffusion coefficients for slab and 2D 
geometry, typical parameter for heliospheric conditions are applied and the remaining wavenumber integrations are 
solved numerically. The numerical calculations are performed for and as given by equations (I82II and (19011 . 
respectively, for the damping model of the DMT approach. 

For the numerical computation, the power spectra 

g{k\i) = C{q)XsiabSBgi^ij{l + (95) 

and 

g{k^) = C{q)X 2 DSBla{l + (96) 

are used for the slab and 2D contribution, respectively. Here, C{q) = {2^/n)~^r{q)/r{q— 1/2), with q = 5/6 being the 
spectral index. For simplicity, the latter is assumed to be equal for slab and 2D geometry. The remaining parameter 
in equations inH) and are the energy densities in slab and 2D fluctuations, and SB^^, respectively, and 

the bend-over scales, denoted by Xsiab and X 2 d- They are proportional to the respective turbulence correlation lengths 
^slab and l 2 D- 

Concerning the correlation lengths, is assumed that I 2 D = ^sZab/10 = 10®m. Furthermore, it is used that the net 
turbulent energy, SB^ = has only a small fraction in slab turbulent energy (20%) and is dominated by 

the 2D contribution (80%), which is believed to be consistent with solar wind observations (Bieber et al. 119941) . The 
background magnetic field Bq is given by 4 • 10“® G, and the ratio 6 B/Bq is choosen to be 0.2. For the decorrelation 
rate Vc, expression (PI is employed (with am = 1 and va = 50 km s ^) for both the slab and the 2D contribution. 
However, in general, different parallel and perpendicular decorrelation processes might govern particle diffusion. 

Numerical results are shown in Figure|21 The dotted line gives the computations for equation II82II . Solutions 
for kY’, Eq. (PTO . are visualized by the dashed and solid curves for electrons (e ) and protons (p), respectively. At a 
glance, for the parameter used here, quasilinear perpendicular diffusion is governed by the transverse structure of the 
turbulent magnetic field. 

However, in view of recent numerical results by Giacaloni & Jokipii (II 99911 (see their Fig. 7 for composite turbulence) 
and Matthaeus et al. cons for a fully three-dimensional and static turbulence, one becomes aware of the fact that 
QLT can not explain their simulations, in contrast to the NLGG theory, where the transverse complexity of the 
turbulent magnetic field plays a crucial role. For static conditions, the NLGG theory provides for a nonvanishing, 
finite diffusion coefficient (see Matthaeus et al. EMI Eq. [7] with 7 (k) = 0). This is not the case in QLT, as it is 
shown in section E and E for slab and 2D turbulence geometry, respectively, since vanishes and 00 for the 

magnetostatic case. This implies for a fully three-dimensional turbulence k± = -I- 00 . QLT can, therefore, 

not provide for an adequate description for perpendicular diffusion, at least for static turbulence. As shown above, 
QLT leads to a nonvanishing, finite diffusion coefficient only if the turbulence is non-static. In QLT, parallel and 
perpendicular diffusion processes are considered as being independent, and the success of the NLGG model is based 
on the nonlinear coupling of these two processes. It becomes clear that such nonlinear effects are crucial for a more 
realistic understanding of spatial particle diffusion and other processes such as particle drift. 
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Another issue concerns the transport of charged particles in turbulent fields having at least one ignorable coordinate. 
It was shown by Jokipii et al. and Jones et al. ifTTina that perpendicular particle diffusion is then suppressed 

and that a particle is tied to the background magnetic field. This result implies that it is not possible to get cross-field 
diffusion for neither slab nor 2D turbulence geometry. The numerical simulations by Qin et al. Il2002bll were performed 
for a “quasi”-slab, magnetostatic turbulence, leading to the insight that perpendicular diffusion is indeed suppressed 
for such a simplified geometry. The particle transport is then subdiffusive and not a Markovian diffusion process. This 
can not be taken into account by neither the NLGC model nor the QLT calculations presented here, since both are 
based on the (Markovian) diffusion approximation. However, though normal diffusion is explicitly assumed by using the 
TGK formula the calculations in section El show that the QLT perpendicular diffusion coefficient approaches zero, 
for both the plasma wave and the dynamical magnetic turbulence viewpoint, if slab fluctuations are static. Giacalone 
& Jokipii iirnMii have shown that cross-field diffusion is also suppressed for a static, two-dimensional turbulence. This 
is in agreement with Jokipii et al. and Jones et al. iirnnHii . In stark contrast to this is the QLT result for a 

two-dimensional turbulence presented in sectional There, it is shown that > oo for static conditions. 

The calculations presented here casts serious doubts on the applicability of QLT for turbulence topologies different 
from slab geometry. This implies that QLT is inapplicable for obtaining conclusions on the nature of interplanetary 
magnetic fluctuations, being in contradiction to the statement made by Shalchi & Schlickeiser II2()()4II . 

Besides being valid only for weak turbulence {6 B/Bq <C 1) and unperturbed particles orbits, the list of limitations 
of QLT can be extended by the result that it fails for non-slab geometries, at least for a static turbulence. Goncerning 
the dynamical behavior of the turbulence. Appendix |S1 describes briefly a modification of the QLT results given 
above, yielding a nonvanishing and finite diffusion coefficient for an “intrinsic” static turbulence. The crude approach 
used there might probably be of more academic interest, since it is based on the speculation that particles alter the 
dynamical behavior of their scattering agent governing their diffusive transport. However, it emphasizes the importance 
of nonlinear contributions for perpendicular diffusion. 

Acknowledgements. The author thanks R. A. Burger, J. Minnie and M. S. Potgieter for valuable conversations. Partial support 
by the South African National Research Foundation (NRF) and the Deutsche Forschungsgemeinschaft (DFG), SFB 591, is 
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Appendix A: Derivation of k 

To derive the quasilinear diffusion coefficient equation 19011 . for 2D turbulence geometry, one can rearrange some 
terms in the corresponding Fokker-Planck coefficient (IHSli, yielding 


p2D _ 27ri?^ 


dk±g(k±)vc 


(1-A^^) .Pn+iiW) + JLiiW) 2 

9 r2 _L „2 + 


E 




(A.l) 


where z = This corresponds to Eq. (54) of Shalchi & Schlickeiser 1200411 . but note that equation (D is derived 
from the general coefficient TO . The two terms including the n-sum can be subjected to a further treatment. Following 
Shalchi & Schlickeiser TOHIi . the formulas 


•7r/2 

Jl_dW) = J d0Jo{2Wcose)cos{20 [n-l]) 


(A.2) 


and 


7r/2 


j 2_^^(TT) = ;^(-l)"+^ J d0Joi2Wcos0)cosi20[n + l]) 

0 

can be used (see Gradshteyn & Rvzhik ITOtiOII . so that 

■k/2 

/ .W,(2WC„S.,C„S(2,„)C»(2.) 
0 

The first term then reads 

oo 

I E 


(A.3) 


(A.4) 


+ JLiiw) 


t/2 
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z sinh(7rz) 


d0Jo{2W cos 6*) cos(26*) cosh(20z) 


(A.5) 
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where the following relation (Gradshteyn & li.vzhik flUtitill is used: 


^ (—1)" /„/, N TJ" cosh(2d2;) 

^4^cos20n = ^ ^ 

A—^ r-* -U 71^ 


1 

2z sinh(7rz) 2z‘^ 

The second term can be treated in a similar manner, resulting in 

r/2 


(A.6) 


E 


n— — C )0 


JliW) 2 

zsinh(7rz) 


d0Jo(21T cos 6*) cosh(2dz) 


(A.7) 


In order to solve the remaining 0-integrations approximately, Shalchi & Schlickeiser II2()()4II used simplified expressions 
for the power spectrum and the turbulence correlation timescale. This is not done here. Instead, the integration with 
respect to fi is first performed. To do so, equations (IA.5II and EH) are first inserted into EIJ. Making use of equation 
l|23), one obtains 


2D AttvRl 


dk±—^ 


9ik±) 


t/2 


sinh(7rz) 


d9 cosh(20z)/^(0, () 


(A.8) 


with 


i 

dti{d,C) = J d^lJoi2C cos0-\/l - - 2(1 - /i^) cos^ 9] 


(A.9) 


where W = k±EL\/l — = CV^ — /J,2 is used. The /r-integration can be solved analytically (Gradshteyn & Ryzhik 

119(1611 to obtain 


ir/2 tt /2 

^iJ.{dyC) = (1 —2cos^0) J c?(/)sin^Jo(2Ccos 0 sin 0)-I-2 cos^ 0 J d^sin ^cos^ 0 Jo(2Ccos 0 sin 0) 


= (1 — 2cos^ 0) 

1 


2n^ r(l/2) (2C cos 0) -b cos^ 0 , Ji (2C cos 0) 


(Ceos 0)3/2 2 


4^3 cos 0 


y/4C COS 9 2 

(1 — 2C^ cos(20)) sin(2C cos0) — 2^ cos0cos(2C cos 9) 


(A.IO) 


where r(a:) and Ji_|_„(a:) denote the Gamma function and spherical Bessel functions of the first kind, respectively. 
With this, the diffusion coefficient (Toil can be written as 


^ 2 d ^ J dk±gikj_)I{C,z) 

0 

with the function z) given by equation (Fmi . 


(A.ll) 


Appendix B: Modification of QLT Calculations 

In view of the failure of QLT in explaining perpendicular diffusion for a static turbulence, the question arises whether 
the quasilinear calculations presented above can be modified to obtain a nonvanishing, finite diffusion coefficient. For 
this, the general perpendicular diffusion coefficient is rewritten, and the DMT viewpoint (damping model) is used. 
Upon substituting the general Fokker-Planck coefficient (ISSll into K ±, Eq. (ISSt . and using the DMT resonance function 
CZHl), one obtains 


K± = 


2 Bl 



Vc _ 

+ (fc||U|| -I- nr2)2 


t/(k,^,n) 


(B.l) 


where Vc = denotes again the turbulence decorrelation rate. Furthermore, 
C/(k, n) = [((fc||UM + nnfv-^ + /i2fc2) Jliw) + (1 - ^^'^)kl [J'^{W)\ 


(B.2) 
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Here, Eqs. iHTl) . lO and m were evaluated for = 0 and Fj = 0 (DMT approach) and inserted into equation 
It was assumed that a — 0. 

Speculated that a number of possible effects may contribute to the decorrelation of magnetic fluctuations. Particles 
diffuse in parallel and perpendicular direction and interact with the turbulence and, therefore, have probably an 
influence on the dynamics of their scattering agent governing their diffusive transport. Proceeding, it is assumed that 
the additional “diffusive feedback” of the particles on the turbulence is given by K||, the parallel diffusion coefficient, 
and by k± itself. This leads for the net turbulence decorrelation rate to the Ansatz 

iyc = '^’^c,i = J + K\\kj + K±k‘i +... (B.3) 

i 

where 7(k) now represents the intrinsic turbulence decorrelation rate. The dots denote other processes such as particle 
drift and stochastic acceleration which might also influence turbulent decorrelation and, therefore, perpendicular 
diffusion. In what follows, equation can be written as 


= 


2 Bl 



_ U{Kk-,n) _ 

(?;/A||)A + K||fcjj + + 7 


(B.4) 


where the following auxiliary function is introduced: 


A(^,n) 


3(^ + n/(fc||i?i))2 

1 + (Kx/K||)(A;_L/fc||)2 +7/«;||A:| 


(B.5) 


It contains the QLT limit of an unperturbed particle orbit. Here, the relation «;|| = A||u/3 has been used, with Ay 
being the mean free path for parallel scattering. Apart from the terms appearing in A, equation IB.411 reveals the same 
nonlinear structure as the recent NLGC result by Matthaeus et al. II2()()3II (see their Eq. [7]). 

Going to the last extrem, it is assumed that particles move only forward and backward to the mean magnetic field. 
This can be taken into account in equation IIB.411 by considering the cases /r = ±1 in terms of a Dirac delta distribution. 
Since W = 0 for ^ = ±1, the Bessel functions in Eq. are nonvanishing only for n = 0. The evaluation then yields 

^ ^ ^ /(k)(l + 2 fcf/fei)/(l + fcf/fci) 

Bq J {v/Xii)A +K\ikf^+K±k'^^+j 


B^J (VA||)+«:||A| + «xfci+7 


(B.6) 


where, for the last step, the denominator is subjected to a rough approximation and A(^ = ±1, n = 0) := 1 is used. A 
closer inspection of Eq. fRAll leads to the result that the diffusion coefficient is nonvanishing and finite for 7 = 0. Apart 
from constants, Eq. IB.611 corresponds to the NLGG result by Matthaeus et al. Il2()()3ll . However, one should keep in 
mind that the derivation of the diffusion coefficient IIB.6II is based on a speculative approach and rough assumptions, 
and that it is presented here for illustrative purposes only. 


References 

Bieber, J. W., & Matthaeus, W. H., 1997, ApJ, 485, 655 

Bieber, J. W., Matthaeus, W. H., Smith, C. W. et al. 1994, ApJ, 420, 294 

Forman, M. A., Jokipii, J. R., & Owens, A. J. 1974, ApJ, 192, 535 

Forman, M. A., 1977, Ap&SS, 1, 83 

Gary, S. P., 2003, Cambridge Univ. Press, New York 

Giacalone, J., & Jokipii, J. R., 1994, ApJ, 430, L137 

Giacalone, J., & Jokipii, J. R., 1999, ApJ, 520, 204 

Gleeson, L. J., 1969, Planet. Space Sci., 17, 31 

Gradshteyn, I.S., & Ryzhik, I.M., 1966, Academic Press, New York 

Jokipii, J. R. 1966, ApJ, 146, 480 

Jokipii, J. R., & Parker, E. N., 1969, ApJ, 160, 735 

Jokipii, J. R., Kota, J., & Giacalone, J., 1993, Geophys. Res. Lett., 20, 1759 
Jones, F.C., 1990, ApJ, 361, 162 

Jones, F.C., Jokipii, J. R., & Baring, M.G., 1998, ApJ, 509, 238 
Kota, J., & Jokipii, J.R., 2000, ApJ, 531, 1067 
Kubo, R., 1957, J. Phys. Soc. Japan, 12, 570 






O. Stawicki: On quasilinear perpendicular diffusion 


19 


Lerche, L, & Schlickeiser, R. 2001, A&A, 378, 279 

Mace, R. L., Matthaeus, W. H., & Bieber,J. W., 2000, ApJ, 538, 192 

Matthaeus, W. H., Goldstein, M. L., & Roberts, D. A., 1990, J. Geophys. Res., 95, 20673 

Matthaeus, W. H., Qin, G., Bieber, J. W., & Zank, G. P., 2003, ApJ, 590, L53 

Qin, G., Matthaeus, W. H., & and Bieber, J. W. 2002a, ApJ, 578, L117 

Qin, G., Matthaeus, W. H., & and Bieber, J.W., 2002b, Geophys. Res. Lett., 29, 1048 

Schlickeiser, R. 2002, Springer-Verlag, Berlin 

Shalchi, A., & Schlickeiser, R., 2004, A&A, 420, 821 

Stawicki, O., Gary, S. P., & Li, H., 2001, J. Geophys. Res., 106, 8273 

Urch, I.H., 1977, AA&SS, 46, 389 



